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Motivated by the successful application of MCRG in momentum space to Acjf> 3 , we determine the critical 
exponents at the crumpling transition in fixed triangulation surfaces. The results are still tentative, but suggest 
that —1.0 > r\ > —1.3, pointing at a value for the fractal Hausdorff dimension at the crumpling transition fixed 
point somewhere between 3 and 4. 



1. INTRODUCTION 

We have to go back to the late eighties (!],§] to 
find the first evidence for the existence of a phase 
transition separating the crumpled and smooth 
regimes in mathematical models of surfaces de- 
fined on a mesh with fixed connectivity. How- 
ever, the actual values of the critical exponents 
associated to this phase transition remain largely 
unknown. Their determination is crucial to iden- 
tify which type of conformal theory one gets in 
the vicinity of the critical point. 

Two recent numerical simulations report the 
following values for the critical exponents a and 
v. In || a direct measuremt of the correlation 
length gives v = 0.71 ± 0.05 and a = 0.58 ± 0.1 
(the latter is obtained via scaling relations). A di- 
rect measurement of the finite size scaling of the 
specific heat gives v = 0.73 ± 0.06, in agreement 
with the previous values. In Q finite size scaling 
gives u = ajv = 1.11±0.1, that is a = 0.71±0.05, 
v = 0.64 ± 0.02. However these values do not 
agree with a direct fit to the specific heat results 
on large lattices (up to 128^1, which tend to give 
much lower values for a |^|J^|. 

Virtually nothing is known of the other criti- 
cal exponents r\ and 7, which cannot be deduced 
from v and a alone using scaling relations. The 
most interesting one is perhaps 77, the anomalous 
dimension of the field Xi that describes the posi- 
tion of the surface. The fractal dimension dn is 
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related to 77 via the relation du — — 4/77. Some- 
times in random surface theory the field anoma- 
lous dimension is introduced through the tangent- 
tangent two-point function. If we denote the cor- 
responding anomalous dimension by 77 it follows 
that 77 = 77 - 2 and d H = 4/(2 - 77). 

In this note we present our preliminary results 
concerning the above critical exponents obtained 
via the Monte Carlo renormalization group. The 
action we shall study is 
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with the sums extending over all pair of neigh- 
bouring sites and triangles (I, J), tik is the 
normal vector to the K-th triangle. Sites and 
triangles live on a two dimensional triangulated 
surface of fixed connectivity embedded in R 3 . 

2. THE METHOD 

Recently we proposed a method of implement- 
ing the renormalization group ideas directly in 
momentum space 0. It works very well in A</>| 
yielding critical exponents at the Wilson fixed 
point with a precision comparable, if not better, 
to any other technique we know of. 

Extremely long autocorrelation times difficult 
the application of local Monte Carlo algorithms 
to the problem at hand. For a 128 2 system 
autocorrelation times can be as large as 10 6 
sweeps [|| .To cure these difficulties a Fourier ac- 
celerated Langevin algorithm was proposed in || 
and used in random surfaces in e.g. ||]. 
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Fourier acceleration works best if the Langevin 
equation is transformed to /c-space. It is there- 
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momentum space directly. The renormalization 
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group transformation that we use is probably the 
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simplest one; we just discard half the Fourier 
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modes in each direction. The interested reader 
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may wish to consult |Q for a more detailed dis- 
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cussion of the method, here we shall dwell on the 
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details specific to random surfaces. 
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In Fourier accelerated Langevin algorithms one 
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uses the freedom to select the two-point function 
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for the gaussian stochastic noise that best serves 
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the purpose of updating the different Fourier 
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modes efficiently. 
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For a free field theory the choice e(j>) = (p 2 + 
to 2 ) -1 , where to is the inverse correlation length 
of the system, works best. This also seems to be 
a good choice for A0|. It has been suggested ||, 
|j| that an effective action that describes many 
features of random surfaces is S e ff — —xAx + 
\xA 2 x and thus it would seem that 
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with to being again the inverse correlation length, 
would be a sensible choice for e(p) (actually, of 
course, the lattice transcription of the above). 
This expectation is not borne out by actual sim- 
ulations, however. It turns out to be nearly im- 
possible to get good convergence of the method 
if one uses the actual inverse correlation length. 
We have found more convenient to select a rel- 
atively large value for the mass, such as to = 1 
which seems to provide a good balance between 
the need to update efficiently the slow modes and 
the need to keep their stochastic excursions in 
check. Most of the results we present below are 
obtained with this value for to. 

A potentially more serious drawback of the 
method is the following. When we update the 
position of the variable on the i-th site we have 
to compute 
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If it just happens that Aj is small, the field vari- 
able will receive a large kick. Because the re- 
sulting configuration is far from the classical tra- 
jectory, it will be promptly brought back to or- 
der by the Langevin algorithm itself, but such 
an unphysical large fluctuation spoils the statis- 
tical samples. It should be emphasized that this 
is a problem associated to the fmiteness of the 
Langevin step and not an intrinsic difficulty of 
our system. One can indeed easily check that 
reducing the Langevin step reduces these off- 
equilibrium excursions, but there is a limit to the 
reduction of the time step as this also increases 
the autocorrelation times. To alleviate this prob- 
lem we modify the action in the following form 



njnj — > njnj exp(— sAt(Aj 1 + A j 1 )) 
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Table 1 illustrates the dependence of the results 
on the parameter s for fixed values of At. Some 
observables can be calculated exactly and so we 
know when the method reproduces the right re- 
sults. We have also analyzed the volume depen- 
dence of s. 

We will discuss next the choice of At and its 
influence on the autocorrelation time r. The 
longest autocorrelation time is always that of the 
gyration radius, but we present results for a fairly 
characteristic one, namely that of the extrinsic 
curvatute Sec = J2i j{^ — ninj). Figure 1 shows 
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Figure 1. Autocorrelation times for Sec at re = 
0.80, s = 10, to = 1 as a function of (At) -1 in 
units of 10 3 (triangles). The squares show the 
dependence on s of the autocorrelation time for 
At = 6 x 10~ 5 , k — 0.80 and to = 1, and corre- 
spond (left to right) to s = 1, 2, 5 and 10. The 
diamonds show the dependence on s of the au- 
tocorrelation time for At — 8 x 10~ 5 , k — 0.80 
and to = 0.125, and correspond (left to right) to 
s = 1, 10, 15, 20. In all cases the size is 16 2 and 
the statistics 5 x 10 5 sweeps. 



the dependence of t on At. (The dependence on 
s for two values of to is also illustrated — there 
is virtually none.). As expected r ~ (At) -1 . In 
view of what we know about the extremely long 
autocorrelation times in local Monte Carlo algo- 
rithms the present autocorrelation times are truly 
impressive. 

One point that practitioners of the Langevin 
algorithm always have to bear in mind is that 
the method is exact only in the At — * limit. 
Indeed the equilibrium distribution is not the one 



corresponding to S but rather to 

Even after tunning s, some systematic errors re- 
main. Fortunately they can be reduced by work- 
ing with a small enough value of At. Yet At 
has to be large enough so that we get reasonably 
short autocorrelation times and manage to sam- 
ple configuration space. After a careful analysis 
we have settled for At = 8 x 10~ 5 as the opti- 
mal value. The renormalization of the parame- 
ters in S is negligible with this value for At. A 
detailed analysis will be provided elsewhere. (Al- 
ternatively, one could work with a second order 
Langevin algorithm, where errors are at least of 
C((At) 2 ). We have tried this, but found no real 
gain.) 

2.1. The Fixed-point Action 

Now we proceed to discuss the essence of the 
method. To determine the critical exponent v we 
need to determine (see e.g. JtJ) the largest eigen- 
value of the matrix T that linearizes the renor- 
malization group transformation in the vicinity of 
the fixed point associated to the crumpling tran- 
sition. In practice truncation in the number of 
operators is required. We have considered a total 
of nine operators. In continuum notation they are 
xAx, xA 2 x and (dxdx) 2 , and those constructed 



with 

TrK = dndx (7) 

TrA 2 = dndn = S EC (8) 

TrA 3 = (d a nd b n)(d b nd a x) (9) 

TrK 4 = (d a nd b n)(d a nd b n) (10) 

Tr(VA) 2 = AnAn (11) 

TrVAT = d a xd a d b n (12) 



and having dimension d < 2. 

We perform a blocking transformation starting 
with our bare action (which consists of just two 
operators), determine T by measuring the appro- 
priate correlators and compute v. In doing that 
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10 6 sweeps, At = ! 


3 x 10" 5 




per site 


big 


small 


Ml 


m 2 = -3.32 


A = 6.0 




0.272 


0.267 




0.161 


0.160 
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R. S. 


k = 0.80 


s = 10 


<^> 

(A 2 ) - (A) 2 
(Sec) 


1.50 
1.49 
1.35 


1.50 
0.86 
1.24 



Table 3 

# of operators Largest eigenvalue 
1st blocking 



1 

3 
5 



2nd blocking 



1 
3 
5 



1.56 
2.09 
2.11 



1.99 
2.22 
2.23 



one needs the field anomalous dimension r/. 77 is 
determined by demanding that at the fixed point 
the renormalizcd and bare actions coincide. It 
is of course unlikely that the fixed point of this 
renormalization group transformation lies just in 
the plane spanned by the two operators of the 
bare action. However we can assess whether we 
are close to it by adjusting 77 and seeing how good 
the agreement between the renormalized and the 
bare observables is. Table 2 illustrates this point 
in a 32 2 — ► 16 2 blocking. Similar results are ob- 
tained in 64 2 32 2 . For comparison some re- 
sults for X(j)f around the Wilson fixed point are 
also shown. We conclude that we are still at some 
distance from the fixed point and that a second 
renormalization is probably necessary. (If we are 
close enough to the critical surface the flow should 
drive us closer to the fixed point.) 

3. THE RESULTS 

It is a very pleasant surprise that the results 
turn out to be extremely sensitive to the value 
of rj used to match the 'big' and 'small' lattices. 
Our results from running on a variety of values 
of k in the vicinity of the crumpling transition 
with large statistics, both in the crumpled and 
the smooth phases and on a variety of systems 
lead us to conclude that i] = —1.05 ± 0.05. from 
the first blocking transformation. 

Our results from the second blocking are still 
very tentative and we have only results for the 
32 2 -> 16 2 ->■ 8 2 blocking, the final lattice being 



undoubtedly too small. Even so, we get a value 
77 ~ —1.3, not too far from the one obtained in the 
first blocking. From these two results we conclude 
that -1.0 > r] > -1.3 and 3.1 < d H < 4.0. 

As for v, Table 3 shows the result of the block- 
ing 32 2 — > 16 2 — > 8 2 . The second blocking yields 
very stable results. Should we use these values 
we would get v = 0.82, but we shall not even 
quote error bars here as these results are still 
very preliminary. Yet, they show considerable 
promise, specially because the computational ef- 
fort required is relatively small. To mention some 
figures, to get the above numbers requires about 
1.5 hours on a Cray YMP or about 5 hours on a 
SGI Power Challenge L. 
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